Draft version February 5, 2008 

Preprint typeset using 1^1^^ style cmulatcapj v. 10/09/06 



BAYESIAN ANALYSIS OF THE LOW-RESOLUTION POLARIZED 3- YEAR WMAP SKY MAPS 

H. K. ERIKSEN^■^■^ Greg Huey^'-^ ^ A. J. Banday", K. M. G6RSKI^•'^■^ J. B. JEWELL■'•^ 
I. J. O'DwYER*-', B. D. Wandelt^-^" 

(Dated: Received - / Accepted -) 
Draft version February 5, 2008 



o 
o 

>>■ 



(N 
O 



> 

rn 

o 
o 



X 



ABSTRACT 

We apply a previously developed Gibbs sampling framework to the foreground corrected 3-yr WMAP 
polarization data and compute the power spectrum and residual foreground template amplitude pos- 
terior distributions. We first analyze the co- added Q- and V-band data, and compare our results to 
the likelihood code published by the WMAP team. We find good agreement, and thus verify the 
numerics and data processing steps of both approaches. However, we also analyze the Q- and V-band 
separately, allowing for non-zero EB cross-correlations and including two individual foreground tem- 
plate amplitudes tracing synchrotron and dust emission. In these analyses, we find tentative evidence 
of systematics: The foreground tracers correlate with each of the Q- and V-band sky maps individu- 
ally, although not with the co-added QV map; there is a noticeable negative EB cross-correlation at 
£ < 16 in the V-band map; and finally, when relaxing the constraints on EB and BB, noticeable dif- 
ferences are observed between the marginalized band powers in the Q- and V-bands. Further studies 
of these features are imperative, given the importance of the low-£ EE spectrum on the optical depth 
of reionization t and the spectral index of scalar perturbations Ug . 

Subject headings: cosmic microwave background — cosmology: observations — methods: numerical 



1. INTRODUCTION 

One of the most remarkable results in the 3-yr data 
release from the Wilkin son Microwave Anisq t ropy Probe 
(WM AP) experiment (|Hinshaw et all l2007t iPage et all 
was the detection of large-scale E-mode polariza- 
tion at millimeter wavelengths. This was interpreted as 
the theoretically predicted signature of reionization, and 
allowed the WMAP team to set new and tighter con- 
straints on the optical depth of reionization r. In turn, 
the well-known degeneracy between r and the spectral 
index of primordial scalar perturbations was broken. 
The final outcome was a claimed de tection of 1 at a 
statistical significance of almost 3cr (jSpergel et al.ll200"7f ) . 

One should bear in mind, however, the great potential 
for systematics effects in both the temperature and po- 
larization measurements. For instance, the precise level 
of power contribution from unresolved point sources af- 
fects Us directl y. An independent a nalys is of this par- 
ticular issue by iHuffenberger et all ()2006[ ) showed that 
the initial point source amplitude quoted by the WMAP 
team was indeed too high, which biased to low values. 
Similarly, on large scales the likelihood approximation 
used by the WMAP team was biased high (jEriksen et al.l 
l2007t l. which also biased low. After these corrections, 
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the statistical significance of rt<; 7^ 1 dropped to ~ 2(t. 

For polarization the situation may be even more seri- 
ous due to the strong sensitivity of r and Hs on the low-.^ 
EE spectrum, combined with the low signal-to-noise ra- 
tio of the WMAP d ata. Systematic eff ects, both from 
the instrument itself (iJarosik et al.ll2007D and from non- 
cosmological foregrounds (jKogut et al.ll2007f ). are much 
more likely to affect the results, and we are also much 
less likely to detect them. It is therefore imperative to 
carefully check both the data and the analysis methods, 
in order to build up confidence in the final cosmological 
results. In this Letter, we start this task by computing 
the low-^ EE, EB, BB and foreground template ampli- 
tude posterior distributions from the WMAP data. 



2. METHOD 

We use a previously introduced Gibbs sampling frame- 
work (Jewell et al. 2004, Wandelt et al. 2004, Eriksen 
et al. 2004, Larson et al. 2007; hereafter JWEL) to esti- 
mate the posterior distributions. For full details on the 
method, we refer the interested reader to the quoted pa- 
pers, and only summarize the principles here. 

First we define our notation. The desired distribution 
is denoted P(s, C^, f |d), where s is the CMB signal, Ce = 
{Cf^, Cf^, Cf^} is the CMB power spectrum, f is a set 
of foreground template amplitudes, and d are the data. 

The Gibbs sampler is a Markov Chain Monte Carlo 
method, and, as such, maps out the full posterior by 
drawing samples from it. While direct evaluation or sam- 
pling from the posterior P{C(\d) requires inversion of a 
prohibitivel y larg e matrix, the Gibbs sampling scheme 
(Gelf and & Smithl ll990) uses the conditional densities of 
the joint posterior P(Cf , s|d) which is computationally 
feasible to sample from. The algorithm may thus be de- 
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Fig. 1. — Comparison of single-^ EE marginal posteriors from the WMAP code (black distributions) and the Gibbs sampler (red 
distributions) using the QV map. The EB and BB power spectra are held at zero. 



scribed by the following sampling steps, 

s'+i^p(s|ci,r,d) 



£4 + 1 



(1) 

-p(Q|s'+i,r,d) (2) 

-P{i\s^+\Cl+\d). (3) 

Here the symbol <— indicates sampling from the condi- 
tional distribution on the right hand side, which can be 
accomplished without inverting the signal-plus-noise co- 
variance matrix (see JWEL for details). For the fore- 
ground template amplitude distribution, we note that 
the required algorithm is identical to that employed for 
sampling monopole and dipole amplitudes (jEriksen et al.l 
[200l . 

3. DATA 

We consider only the low-resolution foreground- 
corrected 3-yr WMAP polarization data in this Letter, 
as provided on LAMBDA^". These come in the form of 
three HEALPix-'^^ sky maps, pixelized at A^sido = 16, each 
having 3072 pixels in both Stoke's Q and U. The WMAP 
P06 sky cut is imposed on the data, leaving only 2267 pix- 
els for the analysis. Two frequency bands are included, 
namely Q-band (41 GHz) and V-band (61 GHz). In ad- 
dition, we analyze the co-added map (denoted QV) , and 
also the two frequency maps jointly but not co-added (de- 
noted Q-fV). All maps are provided with a full noise co- 
variance matrix (jJarosik et al.l [20071) . appropriately cor- 
rected for the P06 sky cut and removal of foreground 

http://lambda.gsfc.nasa.gov 
http://healpix.jpl.nasa.gov 



templates. The units used in this paper are thermody- 
namic fiK. 

For foreground marginalization, we adopt two individ- 
ual templates. First, we use the K-Ka difference map, 
smoothed to 10° FWHM resolution to reduce noise con- 
tributions, as a tracer of synchrotron emission. Second, 
for dust emission we adopt the low-noise template devel- 
oped by the W MAP team for t heir foreground correc- 
tion procedure (iPage et al.ll2007t ). Note that the specific 
shape of these templates are of minor importance; if the 
provided sky maps are free of foregrounds, they should 
not correlate significantly with any non-CMB map. 

We compare our results to the official WMAP likeli- 
hood code^^, also available from LAMBDA. To map out 
the appropriate posteriors, we have written a special- 
purpose MCMC wrapper around this likelihood code. 

4. RESULTS 

4.1. Numerical verification of the WMAP likelihood 

The first case considered is that adopted by the 
WMAP likelihood code, namely the co-added QV map. 
For this analysis, we fix the EB and BB spectra to zero, 
and map out the corresponding marginalized EE poste- 
riors £-hy-£, both with the Gibbs sampler and by the 
WMAP-based MCMC code. 

The results from this exercise are shown in Figure [TJ 
The agreement between the two approaches is very good, 
and this is an important validation of the WMAP data 
processing method: First, we analyze the data at their 

^■^ version v2p2p2 
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Fig. 2. — Marginal posterior distributions for the EE bin power 
between I = 2—6. 



TABLE 1 

Marginalized EE and BB band powers 
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Note. — Values indicate either the posterior mode and 
upper and lower 68% confidence interval, or the upper 68% 
confidence limits. If the lower error bar equals the posterior 
mode value, a non-zero peak is detected, but at a significance 
of less than 68%. 



native A'side = 16 resolution, while the WMAP team 
analyze maps downgraded to TVsido = 8. Second, they 
marginalize over a single total foreground template, while 
we marginalize over the K-Ka difference map and a dust 
template. Third, we use a Gibbs sampler for the nu- 
merical work, while the WMAP team uses a brute-force 
likelihood evaluator. None of these differences affects 
the low--^ EE spectrum peak visibly, as will be quantified 
more precisely in the next section. 

4.2. Generalized analysis 



Having validated the numerics and data processing 
steps in both procedures, we now expand the analysis and 
allow for non-zero EB and/or BB spectra. We also com- 
pute the posteriors for each frequency band separately 
and jointly. For presentational reasons, we report only 
band powers in EE and BB between £ = 2-6 and i = 2- 
20, respectively. Cases of special interest are treated sep- 
arately in subsequent sections. 

In order to achieve good convergence, 10^ samples were 
generated for the non-zero EB, £ = 2-20 cases. For all 
other cases, 10^ samples were generated. The CPU time 
to generate one sample was ~ 2 seconds. 

The results from these computations are summarized 
in Table [T] The EE posteriors with fixed EB = BB = 
are shown in Figure [2l Again, note the excellent agree- 
ment between the WMAP results and the QV and Q-l-V 
cases in the two top sections. 

However, even though the joint QV analyses agree well, 
the picture is considerably less clear when it comes to 
single bands and relaxed EB constraints. On the one 
hand, there appears to be more EE power in the V-band 
data than in the Q-band data. On the other hand, there 
may seem to be small traces of BB power in the Q-band 
data. 

The BB posteriors develop a peak away from 
in the EB^^ case. This is not surprising. Since the 
signal covariance matrix is positive definite, one must 
have C™C^^ > (C™)^. Therefore, when marginahzing 
over C™, a non-zero peak emerges in both the EE and 
BB spectra individually. 

4.2.1. Foreground amplitude posteriors 

In Figure[3]we show the foreground template amplitude 
posteriors for the Q-, V- and QV-band data, for the case 
with fixed EB — BB = 0. Although no signs of significant 
residual foregrounds are observed in the co-added QV 
band, apparently confirming the fits made by the WMAP 
team, the same is not true for each band individually. On 
the contrary, non-zero correlations are seen in both the 
Q- and V-band data individually. 

For the Q-band data, the marginal best-fit K-Ka am- 
phtude is As = -0.027to;"7**, different from zero at 
2(7. The best-fit dust amplitude is Ad = 
For the V-band data, the best-fit dust amplitude is 
Ad = -24.lti'j'i> 2-3cr away from zero. The K-Ka am- 
plitude is As = Omitoml- 

These results may be compared to Table 4 of 
iPage et al.l (|2007D . The main difference between the two 
analyses is that while we apply the conservative P06 
mask to the data. lPage et al.l (j2007l ) apply t he much more 
aggressive processing mask described by iJarosik et al.l 
pOOTi) . The two analyses remove 26.4% and 5.7% of 
the sky, respectively. Considering that all cosmologi- 
cal analyses are carried out with the P06 mask, and 
that variations in the synchrotron spectral index are ob- 
served between the galactic plane and the high latitudes 
(jKogut et al.ll2007h . it is not immediately clear to us why 
the more aggressive mask was chosen for this task by the 
WMAP team. The improvement in raw after further 
correcting the "cleaned" WMAP maps for these residuals 
is -5.4 for the Q-band and -3.3 for the V-band, respec- 
tively. 

Two other, minor differences are that we take into ac- 
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Fig. 3. — The foreground template amplitude marginal posteriors 
spectra were set to zero. 
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Fig. 4. — The ExB cross-power spectrum. The gray region shows 
the la confidence region around the V-band spectrum (blue curve); 
note that the marginal EB posterior is strongly non-Gaussian with 
heavy tails. The green data points show the WMAP EB spectrum 
computed with a pseudo-C^ estimator with Gaussian error bars. 

count both the fuU noise and CMB covariance matrices, 
while the WMAP team includes only a white noise co- 
variance approximation. Note that this carries no extra 
cost within the Gibbs sampling framework. 

4.2.2. ExB cross-correlation spectra 

We find evidence for a non-zero EB correlation in the 
V-band sky map. This spectrum is shown in Figure[4]for 
the Q- and V-band sky maps individually, together with 
the pseudo-C^ EB QV-band spectrum computed by the 
WMAP team. Note the consistently negative correlation 
observed in the V-band. Although the significance of the 
negative correlation is not more than a few tenths per 
multipole, and the joint significance is not more than 1 
to 2(7 depending on binning scheme, it is observed consis- 
tently in every multipole up to ^ = 17. A similar, though 
weaker, trend is seen in the pseudo-spectrum computed 
by the WMAP team from the QV combination. 

5. CONCLUSIONS 

We have performed a Bayesian analysis of the low- 
resolution foreground-corrected 3-yr WMAP polariza- 
tion data using a previously described methodology 




Dust template amplitude 

r the Q-, V- and QV-bands. For this plot, the EB and BB power 

based on the Gibbs sampler. By doing so, we validated 
the numerical implementation of the official WMAP like- 
lihood code, as well as their procedure for degradation of 
map resolution. However, when generalizing the analysis 
to allow for non-zero EB and BB power spectrum compo- 
nents, and also considering individual frequency bands, 
we found several issues that may be important for the 
cosmological interpretation of these data. 

First and foremost, when relaxing the constraints on 
EB and BB, noticeable differences between the Q- and 
V-band posteriors are observed. Specifically, we find gen- 
erally more EE power in the V-band data than in the Q- 
band data, but also perhaps some hints of BB power in 
the Q-band data. At the same time, we have also found 
a negative EB correlation in the V-band map, as well as 
residual foregrounds in both maps. 

If these tentative findings are confirmed by future ex- 
periments or additional years of WMAP observations, 
significant shifts in cosmological parameters could be the 
result. For example, if the V-band data alone were used 
for the WMAP3 analysis instead of the QV combina- 
tion, the amplitude of the £ — 2-6 EE detection would 
increase by 20-50%, depending on whether BB is allowed 
to vary or not. Consequently, r could increase from 0.09 
to ~ 0.12, and rig by a percent or two, comparable to 
its current nominal statistical uncertainty of 0.015, from 
- 0.97 to ~ 0.98. 

As discussed in the introduction, systematics are a se- 
rious concern for both the temperature and polarization 
data for both r and n^. It is important to bear in mind 
that the currently quoted uncertainties on these quanti- 
ties often found in the literature are statistical only. The 
unknown systematic uncertainties may turn out to be 
non-negligible for the currently available data sets, and, 
in particular, we believe it is too early to draw any firm 
conclusions concerning the precise value of n^. Fortu- 
nately, Planck will clarify these issues in the near future. 
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